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DMAHTC  GPS  POINT  POSITIONING  SOFTWARE: 
INITIAL  RESULTS 


THEODORE  E.  MEYER 
GEORGE  T.  TENNIS 
JAMES  A.  SLATER 
BARBARA  J.  DeNOYER 
DEFENSE  MAPPING  AGENCY 
WASHINGTON,  DC  20305-3000  USA 

SUMMARY 

)  Absolute  positioning  software  written  at  the  Defense  Mapping  Agency 
Hydrographic/Topographic  Center  (DMAHTC)  was  designed  to  process  Texas  Instruments 
4100  (TI  4100)  geodetic  receiver  phase  and  range  observations.  This  paper  describes 
absolute  positioning  solutionsobtained  by  processing  data  collected  during  the  Spring  1985 
High  Precision  Baseline  Test  [l]T 

The  positioning  software  can  accept  phase  or  range  observations  and  treat  them  as  either 
pseudoranges  or  range  differences,  thus  allowing  four  different  processing  modes.  The 
software  estimates  receiver  position,  scaling  for  the  tropospheric  correction,  clock  bias,  drift, 
and  aging  parameters,  and  biases  to  the  scaled  phase  measurement. 

Data  from  the  Spring  Baseline  Test  have  been  reduced  and  position  offsets  relative  to 
given  reference  positions  have  been  tabulated.  Comparisons  have  been  made  between 
solutions  computed  using  different  ephemerides,  processing  methods,  observation  types, 
clock  treatments,  and  initial  station  offsets.  Means  and  standard  deviations  of  the  differences 
between  the  position  solutions  and  the  reference  positions  are  presented.  Results  show  that 
the  positoning  software  is  computationally  stable  but  that  significant  day  to  day  variations  of 
the  solution  are  present  due  to  the  data.  (  ' 

1.0  INTRODUCTION 

Absolute  positioning  software  written  at  DMAHTC  was  designed  to  process  TI  4100 
receiver  phase  and  range  observations.  Section  2  of  this  paper  is  a  discussion  of  the 
processing  software  and  methods  used  at  DMAHTC  to  produce  absolute  positions  using 
GPS  data  collected  with  the  TI  4100  receiver. 

Section  3  presents  the  results  of  processing  13  data  sets  collected  at  two  sites  during  the 
Spring  1985  High  Precision  Baseline  Test.  Various  processing  options  were  used  to 
produce  solutions  for  these  data  sets.  Results  are  tabulated  and  compared  to  demonstrate  the 
consistency  of  the  different  processing  options.  Recoveries  of  station  postion  are  shown 
when  the  initial  position  is  offset  by  varying  amounts  in  latitude,  longitude,  and  height. 
Several  data  sets  are  edited  to  restrict  the  satellite  geometry  and  the  results  of  processing  with 
these  data  sets  are  presented. 
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A  flow  chart  of  the  GPS  absolute  point  positioning  process  at  DMA  is  shown  in  Figure 
1.  The  TI 4100  receiver  records  raw  data  collected  from  up  to  four  satellites  at  a  time  onto  a 
cassette  tape.  The  data  on  the  cassette  are  transferred  to  a  9-track  tape  in  exchange  tape 
format  [2]  which  is  then  read  into  a  mass  storage  file  on  a  mainframe  computer.  The 
STARPREP  (Preprocessor  and  Correction)  program  reads  this  mass  storage  file  and  creates 
a  number  of  files  including  the  Corrected  Measurement  File  and  the  Satellite  Transaction  File. 
These  files  and  ephemeris  data  are  combined  in  the  PTPREP  program  to  create  the  input  file 
for  the  GPEST  program.  GPEST  solves  for  the  absolute  point  position  among  other 
parameters. 


Figure  1  -  GPS  Absolute  Positioning  Data  Flow 


All  the  code  is  written  in  modular,  structured  form  using  ANSI  FORTRAN  77  and 
includes  error  recovery  and  execution  trace  features,  restart  capability  for  each  major  section, 
optional  computation  and  application  of  each  of  the  measurement  corrections  in  STARPREP, 
optional  output  of  selected  data,  and  diagnostic  plots. 

2.1.1  SXAREREP  PREPROCESSOR 

The  STARPREP  preprocessor  was  designed  specifically  to  process  GPS  tracking  data 
obtained  from  the  TI  4100  receiver.  The  current  version  of  the  program  requires  input  data 
in  9-track  exchange  tape  format  and  thus  assumes  that  the  original  cassette  tapes  from  the  TI 
4100  have  already  been  converted  to  exchange  tape  format.  It  is  also  implicitly  assumed  that 
the  GESAR  navigation  processor  [3]  was  used  in  the  receiver  to  collect  the  data. 

STARPREP  is  structured  to  handle  one  data  set  (one  tracking  session  for  one  receiver)  at 
a  time  and  leaves  the  function  of  matching  or  combining  related  data  sets  from  multiple 
receivers  to  applications  programs.  There  are  three  major  sections  in  the  program  -  edit,  time 
tag  corrections,  and  data  corrections  -  which  are  executed  sequentially.  The  net  result  is  a 


data  file  containing  measurements  corrected  for  known  sources  of  error  and  the  associated 
corrected  time  tags.  All  of  the  computations  are  done  in  an  earth-centered  earth-fixed  (ECEF) 
reference  system  (WGS-72). 

2.1.1. 1  EDIT  SECTION 

Tracking  data  and  related  information  are  read  into  the  computer  from  the  exchange  tape 
in  the  edit  section  of  STARPREP.  Two  types  of  measurements  are  processed  -  pseudorange 
and  biased  accumulated  Doppler  phase  for  both  LI  and  L2  carrier  frequencies.  Pseudorange 
measurements  are  assumed  to  be  based  on  the  P-code  and  are  scaled  from  the  time  units  of 
the  original  measurement  to  distance  units  by  multiplying  the  uncorrected  pseudoranges  by 
the  speed  of  light.  The  phase  measurement  is  the  continuous  accumulated  cycle  count 
recorded  by  the  receiver,  starting  from  the  time  the  receiver  successfully  "locks  on"  to  a 
given  satellite.  The  phase  measurements  include  an  initial  cycle  ambiguity  (bias)  representing 
the  existing  receiver  cycle  count  at  the  moment  the  local  receiver  clock  is  synchronized  with 
GPS  time.  The  phase  measurements  must  be  corrected  for  an  accumulating  count  due  to  the 
fixed  frequency  bias  purposely  incorporated  in  the  receiver.  To  make  this  correction,  the 
nominal  frequency  offset  is  multiplied  by  the  elapsed  time  from  the  start  of  the  first  good 
observation  after  receiver  synchonization  with  GPS  time.  After  applying  this  term  to  each 
phase  measurement  and  multiplying  the  resulting  quantity  by  the  ratio  of  the  speed  of  light  to 
the  carrier  frequency,  the  phase  measurement  is  scaled  from  units  of  cycles  to  units  of 
distance.  This  "biased  Doppler  range"  value  is  then  used  for  the  rest  of  the  preprocessing. 

The  edit  section  checks  the  data  for  loss  of  signal  or  loss  of  phase  lock,  carrier 
signal-to-noise  ratio  less  than  input  threshold  value,  missing  LI  or  L2  measurements  (time 
gaps),  gross  discontinuities  in  range  or  phase,  and  verification  of  station  position  and 
receiver  equipment  used. 

LI  and  L2  measurements  are  combined  to  compute  two-frequency  ionospheric 
corrections  for  both  the  pseudorange  and  biased  Doppler  range  data.  Only  the  corrected  LI 
data  are  processed  in  the  rest  of  STARPREP. 

In  addition  to  the  measurement  data,  the  edit  section  saves  meteorological  data,  station 
information,  receiver  equipment  information,  and  the  satellite  navigation  message  ephemeris 
and  clock  correction  terms  for  later  use  by  the  preprocessor. 

2. 1.1. 2  TIME  TAG  CORRECTION  SECTION 

The  uncorrected  time  tag  to  each  measurement  is  the  receiver  time  at  the  time  of 
reception.  This  time  tag  must  be  corrected  to  the  time  of  transmission.  This  is  accomplished 
by  dividing  the  measured  pseudorange  by  the  velocity  of  light  and  subtracting  this  value  from 
the  reception  time.  (Propagation  delay  and  the  receiver  clock  offset  from  the  satellite  clock  are 
accounted  for  by  this  simple  process.) 
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An  additional  time  tag  correction  to  account  for  the  satellite  clock  offset  from  GPS  time  is 
made  using  the  navigation  message  values  for  the  satellite  clock  parameters.  This  references 
the  measurements  to  very  nearly  GPS  time. 

2.1. 1.3  DATA  CORRECTIONS  SECTION 

In  addition  to  the  ionospheric  correction  that  was  computed  in  the  edit  section,  five  other 
corrections  are  computed  and  applied  to  the  measurements.  These  corrections  account  for 
tropospheric  delay,  relativity,  satellite  clock,  earth  rotation,  and  satellite  antenna  offset 
The  Hopfield  model  [4]  was  used  for  the  tropospheric  correction  in  conjunction  with  weather 
data  recorded  at  the  site.  Periodic  relativistic  effects  [5]  and  the  offset  of  the  center  of  mass 
of  the  satellite  from  the  satellite  antenna  electrical  center  are  accounted  for.  Measurements  are 
adjusted  for  satellite  clock  offset  and  drift  using  parameters  obtained  from  the  navigation 
message.  Funher  residual  corrections  are  estimated  in  GPEST.  Since  STARPREP  uses  an 
ECEF  reference  system,  an  earth  rotation  correction  is  applied  to  the  data  as  well.  No 
correction  is  made  for  the  receiver  clock  in  the  preprocessor.  For  computing  these 
corrections,  either  the  navigation  message  ephemeris  or  the  Naval  Surface  Weapons  Center’s 
(NSWC)  reference  ephemeris  [6]  can  be  used  as  a  source  of  satellite  positions. 

2. 1.1.4  STARPREP  (TJTPIJT 

The  Satellite  Transaction  File  contains  the  sets  of  orbital  elements  that  are  included  in  the 
navigation  message  for  the  satellites  observed  during  the  time  period  of  the  data  collection. 

The  Corrected  Measurement  File  contains  the  corrected  observations  and  the  value  of 
each  of  the  corrections  as  described  in  Section  2. 1.1.1. 

2.1.2  PTPREP 

The  PTPREP  (Point  File  Preparation)  Program  uses  the  Corrected  Measurement  File  and 
the  Satellite  Transaction  File  to  create  a  Point  File  which  is  used  as  the  primary  input  to 
GPEST.  The  Point  File  consists  of  a  series  of  records  containing  the  observation  time,  the 
corrected  observation,  the  corrections  that  were  applied,  and  a  satellite  position.  The  satellite 
position  may  be  computed  using  the  navigation  message  elements  from  the  Satellite 
Transaction  File  or  interpolated  from  the  NSWC  reference  ephemeris. 

2.1.3  GPEST  ABSOLUTE  POSITIONING 

GPEST  is  an  interactive  absolute  positioning  program  designed  to  accept  data  Files  which 
have  been  processed  by  the  STARPREP  and  PTPREP  programs.  GPEST  consists  of  a 
driver  (GPEST/B1)  and  41  subroutines.  The  program  size  is  about  60k  words  but  can  be 
expanded  if  more  parameters  are  used.The  program  can  use  either  the  pseudorange  (PR)  or 
the  accumulated  Doppler  phase  (AD)  data  type  collected  from  the  TI  4100  receiver.  These 
can  be  processed  using  two  methods;  range  (R),  or  range  difference  (RD).  This  option 
allows  GPEST  to  run  in  four  modes:  pseudorange  processed  as  range  (PR/R),  pseudorange 


processed  as  range  difference  (PR/RD),  accumulated  Doppler  processed  as  range  (AD/R), 
and  accumulated  Doppler  processed  as  range  difference  (AD/RD). 

The  primary  outputs  are  the  adjusted  tracking  station  coordinates  in  WGS-72.  Other 
outputs  include  corrections  to  the  satellite  clock  parameters,  receiver  clock  parameters,  a  scale 
factor  for  the  tropospheric  refraction  correction,  and  the  bias  to  the  phase  count  for  each  span 
of  data. 

Clock  terms  may  be  determined  in  a  global  (gl)  mode  or  in  a  satellite-specific  (sv)  mode. 
"Global"  means  that  only  one  set  of  clock  parameters  (bias,  drift,  and  aging)  occur  in  the 
state  vector  and  that  the  terms  represent  the  receiver  clock  behavior  with  respect  to  GPS  time. 
This  assumes  that  the  satellite  clock  parameters  in  the  navigation  message  correct  each 
satellite  clock  exactly  to  GPS  time.  "Satellite-specific"  means  that  a  set  of  clock  parameters 
appear  in  the  state  vector  for  each  satellite.  A  satellite-specific  term  will  represent  the  relation 
between  the  receiver  clock  and  the  preprocessor  corrected  satellite  clock  . 

To  account  for  biases  introduced  into  the  phase  data  due  to  the  initial  phase  ambiguity  in 
the  TI  4100  or  losses  of  lock,  GPEST  divides  each  data  set  into  "spans."  A  new  span  is 
defined  after  each  discontinuity  (loss  of  lock  or  cycle  slip).  GPEST  uses  the  PR 
measurements  as  a  check  for  large  changes  in  bias  in  the  AD  measurements.  In  the  AD/R 
mode  an  intial  difference  is  computed  between  the  PR  and  AD  measurements  for  a  satellite. 
All  subsequent  measurements  are  likewise  differenced  and  compared  to  the  intitial  difference. 
If  the  absolute  value  of  the  difference  of  these  values  exceeds  a  maximum  value  input  by  the 
user,  then  GPEST  assumes  that  a  phase  discontinuity  has  occurred  and  a  new  "span"  is 
defined.  A  new  range-bias  term  is  then  added  to  the  state  vector.  The  range  bias  terms  are 
not  needed  in  the  AD/RD,  PR/R,  or  PR/RD  mode. 

The  tropospheric  refraction  correction,  which  is  obtained  from  the  Point  File  for  each 
observation,  is  multiplied  by  a  scale  factor  (a)  and  then  added  to  the  corrected  observation 
supplied  by  STARPREP.  Thus,  if  the  correction  is  adequate,  and  other  error  sources  are 
adequately  modeled,  this  factor  should  be  close  to  zero. 

GPEST  utilizes  a  simple  batch  least  squares  fit.  The  state  vector  is  given  by: 

X  =  (  Xsta,  YSIa,  Zsu,  a,  biassv,  driftsv,  agingsv,  range-biassv  sp  )  (1) 

or  by 

£  =(  Xsta’  Ysia’  ZSU’  a-  biasgi’  driftgi.  agingg,,  range-biassv  sp  )  (2) 

depending  upon  the  clock  model  (gl  or  sv)  used.  Any  of  these  terms  are  optional. 

For  equations  1  and  2: 


X$u,  YJta,  ZsU .  Station  position  in  a  north,  east,  up  (NEU)  system 

a . Scale  factor  for  the  tropospheric  refraction  correction 


biasjv,driftsv,agingsv...  Clock  parameters  for  receiver/satellite  clock  combination 


biasgl,driftgl,aginggl....  Clock  parameters  for  receiver  time  with  respect  to  GPS  time 

range-biassv  sp .  Phase  offsets  for  each  satellite  (sv)  and  span  (sp) 

The  observation-state  relation  is  given  by  the  equation: 

Obs  =  Range  +  (a)  Trop  +  c(bias  +  drift  (At)  +  aging(At)2)  +  (range-bias)  (3) 

where: 

Range  =  [  (xsv-  x)2  +  (ysv-  y)2  +  (z$v-  z)2  ] 1/2  (4) 

and 

x,  y,  z . Station  position  in  an  ECEF  system 

xsv,  ysv,  zsv .  Satellite  position  (ECEF) 

At . Time  from  data  set  epoch 

c . Speed  of  light 

Trop . Tropospheric  refraction  correction 

When  the  RD  mode  of  processing  is  used,  both  sides  of  equation  3  are  differenced 
between  two  consecutive  observations  for  the  same  satellite. 

3.0  OVERVIEW  OF  PROCESSING  RESULTS 

The  data  sets  from  two  stations  included  in  the  Spring  1985  High  Precision  Baseline 
Test  will  be  discussed  in  this  paper.  The  Ft.  Davis,  Texas  and  Richmond,  Florida  data  sets 
comprised  the  most  complete  sequences  available  to  the  authors  at  the  time  of  processing. 
The  data  were  collected  on  days  89  through  95,  1985  (Saturday,  30  March  1985  through 
Friday,  5  April  1985).  Figures  2  and  3  show  elevation  angle-azimuth  plots  for  the  data 
collection  interval  on  day  94  at  the  Richmond  site  and  day  93  at  the  Ft.  Davis  site.  These 
plots  illustrate  typical  geometry  for  most  of  the  data  sets.  Data  from  day  93  were  not 
available  for  the  Richmond  site.  Each  data  set  comprises  from  Five  to  nine  hours  of  tracking 
or  from  1900  to  3000  observations  at  30-second  intervals.  Almost  all  data  sets  include 
observations  from  six  satellites.  Day  92  at  the  Ft.  Davis  site  includes  observations  from  only 
five  satellites.  This  data  set  consists  of  only  1319  observations  and  does  not  span  the  same 
daily  time  period  as  the  other  data  sets. 

For  the  processed  data  sets,  PR  measurements  were  generally  on  the  order  of  from 
20000  to  25000  kilometers.  AD  measurements  scaled  in  STARPREP  consisted  of  values 
that  ranged  from  50000  to  -1470000  kilometers.  Using  the  RD  processing  method,  these 
observations  ranged  from  0.01  to  58  kilometers  per  30-second  observation  interval.  The 
noise  on  these  measurements,  which  is  shown  in  residual  plots  produced  by  GPEST,  ranges 
from  1  to  10  meters  for  PR  observations  and  from  5  to  15  centimeters  for  the  AD 
observations.  Thus,  noise  to  measurement  ratios  varied  in  GPEST  from  approximately  10'9 
for  the  AD/RD  mode  to  10"4  for  the  PR/RD  mode. 

The  positions  used  as  initial  estimates  for  computing  corrections  will  be  referred  to  as 
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reference  positions.  The  WGS-72  positions  used  as  initial  estimates  were  transformed  from 
published  NAD  values  [7]  and  geodetic  ties  to  the  GPS  sites.  Richmond’s  initial  WGS-72 
coordinates  were  25T6140657N  latitude,  279T615796E  longitude,  and  -0.022663 
kilometers  in  height.  Ft.  Davis'  initial  WGS-72  coordinates  were  30\6359547N  latitude, 
256*.052905E  longitude,  and  1.58238  kilometers  in  height. 

Tables  1  through  6  and  8  through  10  show  the  approximate  offsets  in  latitude,  longitude, 
height,  and  radial  distance  (in  meters)  between  the  reference  positions  and  the  positions 
computed  by  GPEST  .  If  the  solution  failed  to  converge  after  12  iterations,  "NC"  is  shown. 

If  data  were  not  available,  the  table  entry  contains  " . The  table  columns  represent 

subsequent  days  at  a  site  and  the  rows  represent  the  four  processing  modes.  Tables  1 
through  6  show  the  offsets  produced  using  various  options  in  processing.  These  options 
include  processing  mode,  ephemeris  type,  and  clock  parameters.  Table  7  shows  means  and 
standard  deviations  computed  for  each  coordinate  in  each  mode  for  the  series  of  days  at  each 
site.  Tables  8  and  9  show  only  the  AD/R  and  AD/RD  modes  and  represent  the  offsets 
produced  when  using  initial  positions  that  have  been  perturbed  in  latitude  and  longitude  or  in 
height.  Table  10  shows  offsets  computed  from  subsets  of  several  of  the  original  data  sets. 


Data  sets  collected  over  7  consecutive  days  were  processed  in  the  four  modes  to  test  the 
repeatability  of  each  of  the  methods.  These  results  are  shown  in  Tables  1  through  6  and 
resulting  statisics  are  presented  in  Table  7. 

Two  measurements  types,  PR  and  AD,  were  collected  at  30-second  intervals  by  the  TI 
4100.  These  data  types  were  described  in  Sections  2. 1.1.1  and  2.1.3.  Both  types  are 
treated  in  much  the  same  manner  by  GPEST.  However,  the  AD  type  requires  the  estimation 
of  an  additional  bias  (see  Section  2.1.3)  for  each  "span"  of  satellite  data.  Tables  1  through  6 
show  the  results  of  processing  using  both  types  of  observations. 


Table  1-  Solution  Offsets  for  Richmond  Site  Using 
NSWC  Reference  Ephemeris  and  Satellite-Specific  Clock 

Meters  from  Reference  Position _ 


Dav  89 

Dav  90 

Dav  91 

Dav  92 

Dav  93 

Dav  94 

Dav  95 

AO 

+  3.59 

+  2.13 

+  3.4  1 

+03.37 

+  2.63 

+  0.57 

AX 

+  5.89 

+  4.4  1 

+  9.82 

+  1  1.70 

+  3.74 

+4.39 

< 

mm 

-5.99 

-2.35 

+  0.63 

-01.93 

+  1.50 

+  5.66 

Rad 

+  9.11 

+  5.44 

+  10.43 

+  12.34 

+  4.82 

+  7.19 

AO 

+  2.10 

+  1.75 

+  0.74 

+  01.55 

+  2.43 

+  2.14 

'**2 

y*. 

+  S.08 

-6.92 

+  8.66 

+  07.99 

+  3.43 

+  4.59 

< 

— 

A  h 

+  0.03 

-0.54 

+  1.61 

-00.28 

-1.56 

-0.62 

Rad 

+  8.37 

+  7.17 

+  8.86 

+  08.17 

+  4.52 

+  5.12 

AO 

+  3.58 

+  3.56 

NC 

NC 

+  2.8  1 

+  0.15 

V 

A  A. 

+  7.65 

+  5.23 

NC 

NC 

+  4.77 

+  6.52 

— 

■til 

-0.62 

+  2.03 

NC 

NC 

+  4.64 

+  5.43 

Rad 

+  8  .as 

+  6.65 

NC 

NC 

. 

+  7.23 

+  8.50 

_ 

AO 

A/u 

+  9.97 

+  9.4  1 

-0.55 

-12.35 

-8.85 

-1.98 

+23.85 
+42.1 1 

. . 

+  7.37 

+  11.10 

+  7.12 
+  29.67 

— 

-✓ 

A  h 

-0.42 

-9.59 

-14.38 

-35.26 

-7.06 

+13.10 

Rad 

+13.69 

-15.65 

+16.98 

+48.34 

+  15.08 

+33.23 

As  described  in  Section  2.1.3,  processing  can  also  be  done  using  two  methods,  R  and 


RD.  The  R  processing  method  treats  both  measurement  types  as  biased  pseudoranges.  The 
RD  method  differences  consecutive  measurements  for  a  satellite,  eliminating  biases  which  are 
common  to  a  satellite.  This  processing  method  eliminates  the  need  to  solve  for  the  individual 
span  biases  and  other  unmodeled  biases.  Tables  1  through  6  compare  the  results  of 
processing  using  both  methods.  Tables  1  and  2  show  comparisons  of  results  from  the  four 
modes  of  processing  for  the  Richmond  and  Ft.  Davis  sites.  The  NSWC  reference 
ephemerides  and  satellite-specific  clock  parameters  were  used  to  produce  these  two  tables. 

Table  1  shows  PR  and  AD  measurements  processed  in  the  R  mode  produce  similar 
results.  For  the  six  sets  of  data  processed  from  the  Richmond  site,  the  AD/R  mode  produced 
mean  offsets  of  2.62  meters  for  the  latitude  (0),  6.66  meters  for  the  longitude  (X),  and  -0.41 
meters  for  the  height  (h)  as  shown  in  Table  7.  This  compares  closely  in  <{>  and  X  to  the  PR/R 
mode  offset  means  computed  from  only  four  solutions  as  shown  in  Table  7.  For  PR/R  mode 
means  of  2.53,  6.04,  and  2. 87  meters  for  6,  X,  and  h,  respectively,  were  produced.  The 
AD  RD  mode  shows  similar  results  with  offset  means  of  1.79,  6.61,  and  -0.85  meters  for  <5, 


X.  and  h,  respectively.  The  R  RD  mode  shows  the  largest  differences  producing  station 
offset  means  of  6.49.  17.11.  and  -2.96  meters  for  6.  X,  and  h,  repectively. 


Table  2  -  Solution  Offsets  for  Ft.  Davis  Site  Using 
NSWC  Reference  Ephemeris  and  Satellite-Specific  Clock 

Meters  From  Reference  Position _  _ 


Dav  89 

Dav  90 

Dav  91 

Day  92 

Dav  93 

Dav  94 

ESSESHl 

AO 

+  3.08 

-0.84 

+  1.11 

-0.03 

-0.43 

+  3.14 

D 

AX 

+  7.70 

+12.93 

+12.15 

+  11.92 

+  12.20 

+22.64 

< 

mm 

+  1.64 

+2.91 

+  0.05 

-0.30 

+  0.77 

+  1.69 

-0.35 

ESI 

+  8.44 

+13.27 

+  12.19 

+  6.11 

+  1  1.93 

+  12.31 

+22.88 

AO 

+  1.32 

+  0.60 

-0.39 

+  0.37 

+  0.71 

+  1.35 

a 

D 

AX 

+  6.55 

+  13.25 

+  12.79 

ii-fs 

+  13.32 

+  12.99 

+13.75 

< 

Cxi 

warn 

+  0.31 

+  1.14 

+  0.94 

+  0.96 

+  0.46 

+  0.69 

+  0.46 

Rad 

+  6.68 

+  13.31 

+12.82 

+  7.72 

+  13.33 

+13.02 

+  13.83 

AO 

NC 

+  1.44 

+  1.66 

+  0.32 

+  1.47 

Cxi 

4X 

NC 

+13.54 

+  19.52 

+  12.32 

+  21.35 

CL 

HS 

NC 

-1.45 

+  3.62 

+  0.21 

+  1.94 

Rad 

+  8.50 

NC 

+13.68 

+  19.92 

+  12.32 

+  21.51 

Cxi 

Q 

AO 

AX 

+  6.48 

+  9.65 

- 1 .04 

+18.15 

+  2.71 
+  2.01 

+  3.60 

-4.25 

+  0.59 

+13.18 

+  6.68 
+31.74 

+  11.11 
+  23.56 

CL 

Cxi 

mm 

-3.97 

+  1 .84 

+  0.26 

-23.42 

+  4.46 

+ 1 .02 

+  1.57 

Rad 

-12.26 

+18.26 

+  3.37 

+24.07 

+  13.93 

+  32.45 

+26.09 

The  Ft.  Davis  site  offsets  are  shown  in  Table  2.  Means  and  standard  deviations  for  the  7 


solutions  are  shown  in  Table  7.  The  Ft.  Davis  site  shows  internal  consistency  similar  to  the 
Richmond  site  among  the  offsets  produced  by  the  AD/R,  AD/RD,  and  PR/R  modes.  The 
PR/RD  mode  again  shows  the  largest  divergence  from  the  other  modes. 

The  AD/RD  mode  produced  the  smallest  standard  deviations  in  0,  X,  and  h  for  both  the 
FT.  Davis  and  the  Richmond  sites  in  all  but  one  case.  The  standard  deviaton  in  latitude 
resulting  from  the  PR'R  mode  for  the  Richmond  site  is  smaller  than  that  for  the  AD/RD 
mode.  However,  only  four  solutions  were  used  to  produce  this  statistic  due  to  the  inability 
of  the  PR/R  mode  to  converge  for  several  data  sets.  In  every  case  the  PR/RD  mode 
produced  the  highest  standard  deviations. 


3.1.2  EPHEMERIS  /NAVIGATION  MESSAGE/NSWC  REFERENCE) 

Processing  was  done  for  the  Ft.  Davis  and  Richmond  sites  using  the  two  ephemeris 
types  available  to  the  STARPREP  and  PTPREP  programs.  The  navigation  message 
ephemeris  was  computed  using  the  algorithm  described  by  Van  Dierendonck  [8].  The 
reference  ephemeris  was  produced  by  NSWC  and  was  evaluated  using  a  Lagrange  8-point 
interpolator.  The  reference  ephemeris  has  a  resolution  of  15  minutes. 

Tables  3  and  4  show  solutions  computed  using  the  navigation  message  ephemeris  for  the 
Richmond  and  Ft.  Davis  sites.  These  may  be  compared  to  Tables  1  and  2  which  show 
solutions  using  the  interpolated  NSWC  reference  ephemeris.  Tables  3  and  4  show  results  of 
processing  with  satellite-specific  clock  parameters  as  do  Tables  1  and  2. 


c  qIax 
<  a: 


Dav  89 


+  3.30 
+9.59 
-8.72 


15 

42 

4: 

33 


+  1.61 
+  13.05 
+  8.68 
+15.76 


+  8.03 
+  11.96 
+4.52 
+  15.08 


Table  3  -  Solution  Offsets  for  Richmond  Site  Using 
Navigation  Message  Ephemeris  and  Satellite-Specific  Clock 

Meters  from  Reference  Position 


Dav  94 


Dav  90 


+2.04 

+9.06 

-2.56 

+9.64 


+  1.44 
+  10.03 
-2.27 
+  10.39 


-1.97 
+  8.01 


-0.88 
+16.28 
-1  1.74 
+  20.10 


Day  91 

Dav  92 

+4.25 

+  0.76 

+  12.67 

+  15.46 

+  0.33 

+2.42 

+  13.37 

+15.68 

+0.52 

-2.06 

+  9.69 

+  7.46 

+  0.52 

+  4.11 

+  9.72 

+  8.77 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

-9.62 

+20.20 

+  0.27 

+  24.31 

-14.78 

-28.64 

+17.62 

+42.62 

Dav  95 


+  4.75 


2.82 

5.90 


+  4.21 
+  7.77 


6.63 

2.83 


+  7.28 
+  16.16 


+  5.47 
+  29.08 
+  13.46 
+  32.52 


Table  7  shows  the  means  and  standard  deviations  for  the  solutions  produced  using  both 
ephemeris  types.  The  latitude  offsets  produced  using  the  navigation  message  ephemeris  for 
both  the  Ft.  Davis  and  Richmond  sites  were  smaller  than  those  produced  using  the  NSWC 
reference  ephemeris  for  all  modes.  The  longitude  offsets  produced  using  the  navigation 
message  ephemeris  for  both  the  Ft.  Davis  and  Richmond  sites  were  larger  for  all  modes. 
The  mean  heights  were  affected  differently  by  the  two  ephemeris  types.  The  Richmond 
mean  heights  were  consistently  higher  using  the  NSWC  reference  ephemeris  and  the  Ft. 
Davis  mean  heights  were  consistently  lower. 

For  the  Richmond  site  the  o,  and  h  standard  deviations  were  smaller  for  the  reference 
ephemeris  in  the  AD/R  and  AD/RD  modes.  Ft.  Davis  shows  higher  standard  deviations  for 
the  reference  ephemeris  for  9  and  A.  in  every  mode  except  PR/R.  In  PR/R  mode  the  standard 
deviaitons  for  <t>  and  X  were  smaller  for  the  reference  ephemeris  and  the  standard  deviations 
for  the  heights  were  smaller  for  the  reference  ephemeris  in  all  modes. 
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Table  4  -  Solution  Offsets  for  Ft.  Davis  Site  Using 
Navigation  Message  Ephemeris  and  Satellite-Specific  Clock 

Meters  from  Reference  Position 


Dav  89 

Dav  90 

Day  91 

Dav  92 

ESK9H 

HTTlff 

Dav  95 

AO 

■Rm 

-01.11 

+  1.17 

-0.53 

-0.05 

-0.95 

+  2.24 

AX 

■til 

+  13.80 

+  12.27 

+  8.21 

+  1 1.68 

+  12.46 

+22.90 

< 

-2.01 

-00.79 

-2.17 

-4.47 

-1.69 

+  0.01 

-2.17 

ESI 

+  8.7  1 

+  13.87 

+  12.53 

+  9.37 

+  11.81 

+  12.50 

+23.15 

AO 

+  0.54 

-00.22 

-0.25 

+  0.94 

+0.59 

+0.42 

r* 

AX 

+  7.44 

+14.23 

+13.29 

+  9.17 

+12.75 

+  13.74 

+  12.80 

< 

zc 

09 

-2.43 

-01.85 

-1.42 

-2.95 

-1.10 

-0.18 

-0.70 

Rad 

+  7.85 

-►14.36 

+  13.38 

+  9.68 

+  12.81 

+13.77 

+  12.85 

AO 

+  3.08 

NC 

+  1.47 

+  1.26 

-0.06 

-0.87 

+0.89 

AX 

+  7.82 

NC 

+13.59 

+  21.77 

+  12.07 

+14.77 

+  21.61 

09 

-1.62 

NC 

-3.82 

-0.37 

-2.34 

+  1.61 

+  0.69 

Esa 

+  8.56 

NC 

+14.20 

+  21.28 

+  12.31 

+  14.90 

+  21.68 

AO 

+  6.12 

-01.80 

+2.30 

+  2.93 

-0.01 

+  6.05 

+9.90 

AX 

+10.69 

-19.63 

+  3.33 

-2.03 

+14.12 

+33.09 

+22.37 

— 

2C 

HI 

-7.35 

-01.41 

-2.03 

-26.92 

+2.55 

-0.1  1 

+0.29 

Rad 

+14.34 

*19.78 

+4.52 

+27.15 

+  14.37 

+33.67 

+  24.48 

3.1.3  CLOCK  TREATMENT  ( GLOB  A  L/S  ATELLITE-SPECIFIO 


GPEST  can  solve  for  clock  parameters  in  the  satellite-specific  mode  and  in  the  global 
mode  as  described  in  Section  2.1.3.  Aging  was  not  estimated  in  any  of  the  data  sets 
discussed.  STARPREP  uses  the  navigation  message  clock  parameters  to  correct  the  satellite 
clocks  before  any  clock  parameters  are  estimated  in  GPEST. 

Table  5  -  Solution  Offsets  for  Richmond  Site  Using 
NSWC  Reference  Ephemeris  and  Global  Clock 

_ _  Meters  from  Reference  Position 


Dav  89 

Dav  90 

Dav  91 

Dav  92 

Dav  93 

Dav  94 

ESKSH 

A0 

+  3.42 

-3.04 

+  2.13 

+  2.81 

+  1.68 

/— > 

AX 

+  6.10 

-6.67 

+  9.49 

*8.17 

+  4.01 

< 

-3.4  1 

-2.46 

+  1.56 

+  2.83 

+  6.05 

-1.19 

nna 

+  7.78 

+  7.74 

+  9.87 

+  9.10 

| 

+  7.45 

+  7.19 

A0 

-2.45 

+  1.76 

Min 

+  2.12 

+  1.96 

AX 

-6.38 

+  8.51 

+  3.09 

+4.19 

< 

oc 

031 

-2.47 

+  0.27 

-0.87 

-1.14 

Rad 

+  7.68 

-7.72 

+  8.71 

+  7.87 

+  3.92 

+4.67 

A0 

-6.07 

NC 

NC 

KHHgHI 

+  1.89 

+  3.92 

AX 

1 

-1.61 

NC 

NC 

+  2.01 

+4.04 

Q9 

1  Bill 

-7.35 

NC 

NC 

+  9.44 

+  3.28 

Rad 

-9.66 

NC 

NC 

+  9.84 

+  14.43 

AO 

+  4.54 

wmggm 

-4.43 

-1.42 

+  9.76 

+  7.51 

c. 

AX 

+  3.79 

Hamaoi 

+12.54 

+  3.03 

+  6.15 

2L 

KIM 

+  6.74 

1 

-10.47 

-1.43 

-4.53 

Rad 

+13.10 

+  16.95 

+  19.10 

10.29 

+  10.77 

Solutions  shown  in  Tables  1  through  4  were  produced  using  the  satellite-specific  mode 


for  the  bias  and  drift  terms  in  processing.  Tables  5  and  6  represent  processing  done  with  the 
bias  and  drift  terms  determined  globally. 

The  solution  offsets  presented  in  Table  5  produced  statistics  that  show  improvement  of 


position  consistency  in  every  mode  except  PR/R  when  compared  to  Table  1.  Table  7  shows 
that  the  global  clock  treatment  produces  lower  standard  deviations  in  all  coordinate  offsets 
and  modes  except  for  the  PR/R  mode  and  the  X  offset  in  the  AD/RD  mode.  In  the  case  of  the 
X  coordinate  in  the  AD/RD,  standard  deviations  are  nearly  equal  (2.13  meters  for  the 
satellite-specific  clock  mode  as  opposed  to  2.12  meters  for  the  global  clock  mode). 

Table  6  -  Solution  Offsets  for  Ft  Davis  Site  Using 
NSWC  Reference  Ephemeris  and  Global  Clock 

Meters  from  Reference  Position 


Dav  91  Day  92  Dav  93  1  Day  94  I  Day  95 


I  A6 

a  a! AX 

<  as 


+0.82 
+  8.05 
+2.87 
+  8.61 


-0.02 
+10.50 
+0.97 
+  10.58 


+  9 

.95 

NC 

NC 

NC 

NC 

-1. 

38 

+  14 

.61 

+  3 

.41 

+  15 

.10 

+  0.44 
+  8.25 
-2.70 
+  8.71 


+  10.06 
-0.26 
+10.11 


+  0.30 
+  7.37 
+  5.55 
+  9.24 


+  0.33 
+  4.28 


+7.70 


+3.81 

-3.13 

-1.18 

+5.07 


+  17.92 


+  0.54 

+  0.39 

+  8.01 

+  1 1.08 

-2.91 

+  0.72 

+  8.56 

+11.15 

+0.67 

+  0.85 

+  10.29 

+9.90 

-0.41 

-0.18 

+10.35 

+9.97 

+  0.36 

-0.25 

-1.76 

+  2.19 

-3.47 

-5.66 

+3.91 

+  6.07 

+  3.07 

-0.69 

+13.58 

+23.97 

.  +2.31 

+4.77 

+14.15 

+24.51 

+  0.98 
+9.51 
-6.54 
+  11.60 


-7.83 
+  8.78 


+  1.99 
+15.70 
+6.53 
+17.16 


Table  7  ■ 


A0  (Meters) 


Sv-Specific  Clk  G1  Clk 


Nv  Eph  Reference  Enh 


Summary  of  Offset  Statistics 


AX(Meters)  Ah(Meters) 


Sv-Specific  Clk  G1  Clk  Sv-Specific  Clk  G1  Clk 


Nv  Ephl  Reference  Eph  Nv  Eph  Reference  Eph 


+8.37  +6.04 

+  15.79  +17.11 


AD/RD 


PR/R 


PR/RD 


AD/R 


AD/R 


PR/R 


PR/RD 


AD/R 


AD/R 


PR/R 


PR/RD 


1 .46 
1.46 
1 .21 
9.91 


+  0.50 
+  0.33 
+  0.96 
+  3.64 


6.89 

6.26 


+  0.91 
+  6.93 


10.14  15.91  4.09  16.01  11.61  8.46 


+  1.32 
+4.30 


+12.76  +12.24 
+  11. 92  +  I  1.46 
+15.27  +14.81 
+14.46  +13.43 


r*i«j 


tl»! 


Three  days  of  data  for  both  sites  were  processed  using  initial  position  vectors  that  were 
offset  from  the  reference  postion  in  <J>  and  X.  These  perturb<*Ton<of  the  initial  position  were 
done  in  order  to  test  GPEST's  ability  to  recover  a  position  given  an  inaccurate  initial 
position.  This  is,  in  one  sense,  a  test  of  the  computational  stability  of  the  positioning 
procedure.  The  offsets  were  on  the  order  of  4  arcseconds  or  approximately  110  to  120 
meters.  The  Richmond  site  is  perturbed  by  approximately  110  meters  in  latitude  and  -110 
meters  in  longitude.  The  Ft.  Davis  site  is  perturbed  by  approximately  -110  meters  in  latitude 
and  1 10  meters  in  longitude.  These  figures  might  be  typical  of  the  accuracy  of  coordinates 
scaled  from  a  map  grid. 

Table  8  presents  offsets  produced  with  these  initially  perturbed  latitudes  and  longitudes. 
These  offsets  can  be  compared  directly  to  results  from  Tables  1  and  2.  Aside  from  the  initial 
station  coordinates,  all  modes  and  input  parameters  were  the  same  for  the  results  presented  in 
both  sets  of  tables.  Tne  differences  in  the  solutions  are  minimal.  As  an  example,  the  6,  X, 
and  h  offsets  for  day  90  from  Table  8  in  the  .AD/RD  mode  are  1.75,  6.94,  and  -0.54  meters. 
Table  1  shows  values  of  1.75,  6.92.  and  -0.54  meters  for  the  coordinates.  The  total 
difference  is  only  2  cm.  Tr.e  other  data  sets  shown  in  Table  8  show  similar  small  differences. 

Table  S  -  Solution  Offsets  w-ith  Initial  Latitude  and  Longitude  Perturbed 

NS'-VC  Referer.ee  Eshetneris  -  SateUite-Soecific  Clock  -  All  Uaiis  arc  Meters 


Dav  90 


R-.chrnor.d  : 
1  Dav  ?T 


Dav  92 


-2.03 

-4.69 

-3.35 

-6.-0 


-3.16 
*9.23 
-0.69 
-9. 7g 


-0.73 


-1.53 

-3.35 


-3.37 

-11.71 

-1.93 

-12.34 


-1.54 

-3.0! 

-0.23 

-3.16 


Ft.  Davis  Site 
~1  093  r 


-0.43 
+12.18 
+  1.69 
-12.30 


-0.96 

-7.73 


Tnree  days  of  data  from  both  sites  were  processed  using  initial  position  vectors  that  were 
offset  in  h.  The  offsets  were  10  meters  for  day  90,  30  meters  for  day  91,  and  780  meters  for 
day  92  at  the  Richmond  station  and  -40  meters  for  day  92,  10  meters  for  day  93,  and  -100 
meters  for  day  94  at  the  Ft.  Davis  site. 

Table  9  presents  offsets  produced  with  these  initially  perturbed  heights.  Again,  these 
offsets  can  be  directly  compared  to  results  from  Tables  1  and  2.  As  with  the  latitude  and 


longitude  offsets,  the  differences  in  the  solutions  are  minimal.  As  an  example,  the  <J>,  X,  and 
h  offsets  for  day  90  from  Table  9  in  the  AD/RD  mode  are  1.75,  6.94,  and  -0.54  meters. 
Table  1  shows  values  of  1.75,  6.92,  and  -0.54  meters  for  the  coordinates  for  the  same  day 
and  processing  mode.  As  with  the  latitude  and  longitude  offsets,  the  total  difference  is  only 
2  cm.  The  other  data  sets  shown  in  Table  9  exhibit  similar  small  differences  when  compared 
to  Tables  1  and  2. 

Table  9  -  Solution  Offsets  with  Initial  Heights  Perturbed 

NSWC  Reference  Ephemeris  -  Satellite-Specific  Clock  -  All  Units  are  Meters 


Q  J  AX 
< 


AO 

c  cl  ax 

<  oc 


Richmond  Site 


Dav  90  Dav  91  I  Dav  92 


+  3.16 
+  9.23 
-0.69 
+  9.78 


+  0.74 
+  8.68 
+  1.61 
+  8.86 


+2.30 
+  4.69 
-3.85 
+  6.40 


-0.54 
+  7.17 


+  3.37 
+  1  1.72 
-1.93 
+  12.35 


1.5 

8.0 

0.2 


+  8.18 


092 


-0.11 
+  6.10 
-0.93 
+  6.11 


Ft.  Davis  Site 


093 


+0.37 

+00.71 

+  13.31 

+  12.98 

+  0.49 

+  0.69 

+13.32 

+13.02 

Two  of  the  original  data  sets  were  edited  to  create  subsets  that  could  be  used  to  test  the 
effects  of  various  factors  on  the  station  position  solutions.  Some  of  these  factors  are  the 
amount  of  data,  the  time  span  of  the  data,  the  number  of  satellites  from  which  observations 
are  used,  and  the  change  in  satellite  geometry.  Figures  2  through  7  illustrate  the  geometry  of 
the  data  used  for  these  sets.  The  outer  ring  in  the  figures  represents  an  elevation  angle  of  0”. 
The  inner  circle  represents  45°,  and  the  crosshair  represents  the  zenith  point. 

Table  10  -  Solution  Offsets  for  Edited  Data  Sets 

NSWC  Reference  Ephemcris  -  Satellite-Specific  Clock  -  AU  Units  are  Meters 


Richmond  Site  -  Dav  94 


North  Quads  South  Quads 


C  Cl  AX 
<  =4 


+  2.91 
+  1.85 
+  0.55 
-3.48 

+  7.45 
+23.92 
+  6.36 
+25.90 

-2.90 

+  6.20 

+  1 .89 

+  1  8.72 

-0.69 

+  4.39 

-3.52 

+20.24 

Ft.  Davis  -  Dav  93 


North  Quads  South  Quads 


-0.65  -11.44 

+16.92  -2.45 

-2.82  +14.00 

+17.21  +18.23 


+0.75  -11.93 

+21.28  +1.96 

+1.30  -9.84 

+21.39  +15.57 


Figure  2-Richmond-Day94  Figure  4-Richmond- Day94  Figure  6~Richmond-Day94 

Unedited  Southern  Quadrants  Northern  Quadrants 
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Unedited 


Figure  5-Ft.Davis-Day93  Figure  7-Ft.Davis-Day93 
Southern  Quadrants  Northern  Quadrants 
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The  data  set  that  was  collected  on  day  94  at  the  Richmond  site  was  edited  twice  to 
produce  two  smaller  data  sets,  one  with  observations  limited  to  those  taken  while  the 
satellites  being  tracked  were  in  the  northern  quadrants,  and  another  with  observations  limited 
to  the  southern  quadrants.  The  data  set  for  day  93  at  the  Ft.  Davis  site  was  edited  to  produce 
two  subsets  in  a  similar  manner.  Figures  2  and  3  represent  the  unedited  data  sets.  Figures  4 
and  5  represent  the  data  sets  limited  to  the  southern  quadrants  and  Figures  6  and  7  represent 
the  sets  limited  to  the  northern  quadrants.  Table  10  shows  the  offsets  produced  by  reducing 
these  data  sets.  These  sets  were  processed  using  satellite-specific  clock  parameters  and  the 
reference  ephemerides  and  may  be  compared  to  results  shown  in  Tables  1  and  2. 

In  the  AD/RD  mode,  the  edited  data  set  which  most  closely  recovers  the  offsets 
computed  from  the  Table  1  solutions  is  the  set  shown  in  Figure  6.  This  set  contains  1051 
observations  as  opposed  to  251,  547,  and  472  observations  for  the  other  three  edited  data 
sets.  The  solutions  from  the  other  three  edited  sets  do  not  fall  within  one  standard  deviation 
of  the  solutions  for  the  unedited  sets. 


From  the  solutions  presented,  a  systematic  offset  was  detected  in  latitude,  longitude,  and 
height  for  each  of  the  two  stations  in  all  modes.  The  source  of  this  offset  is  unknown  but 
may  stem  from  several  sources  such  as  error  in  the  reference  positions,  ephemeris  error,  or 


systematic  error  due  to  satellite  geometry. 

There  was  a  significant  amount  of  day-to-day  variation  in  all  coordinates  and  modes. 
However,  as  was  shown  in  Table  7,  the  accumulated  Doppler  count  observation  type 
processed  as  range-differences,  with  global  clock  parameter  solution  and  the  reference 
ephemeris,  produced  standard  deviations  in  <{>,  X,  and  h  which  in  the  majority  of  cases  were 
smaller  than  the  other  modes.  The  offsets  for  this  mode  show  standard  deviations  over  a 
week  which  were  sub-meter  in  latitude  and  height  and  approximately  2  meters  in  longitude. 
The  Doppler  data  processed  as  ranges  and  the  pseudorange  data  processed  as  ranges 
generally  produced  standard  deviations  that  were  only  slightly  less  consistent  then  the 
AD/RD  mode.  The  pseudorange  data  processed  as  range  differences  produced  the  highest 
standard  deviations  for  station  position  in  every  mode.  This  may  be  due  to  the  measurement 
to  noise  ratio  as  described  in  Section  3.0.  For  the  two  sites  longitude  consistently  shows  the 
largest  offsets  in  all  modes. 

The  position  estimation  procedure  in  GPEST  appears  to  be  insensitive  to  perturbation  in 
the  initial  estimates  of  the  station  coordinates  of  up  to  780  meters.  The  results  shown  in 
Tables  8  and  9  demonstrate  that  the  solution  converges  to  the  same  position  regardless  of  the 
initial  offset  of  the  a  priori  coordinates.  More  tests  of  this  nature  need  to  be  conducted  to 
determine  the  limit  of  this  stability. 

Satellite  constellation  configuration  and  the  amount  of  data  obviously  affect  the 
positioning  solution.  Editing  the  data  sets  to  bias  the  geometry  of  the  data  processed  by 
GPEST  has  an  affect  on  the  solution,  but  the  few  cases  explored  here  are  inconclusive  as  to 
any  systematic  effect.  In  the  cases  represented,  the  amount  of  data  and  the  data  interval  may 
have  been  the  determining  factor  rather  than  geometry.  With  short  spans  of  data,  periodic 
ephemeris  errors,  which  are  averaged  out  over  longer  periods,  may  become  more  dominant 

Further  experimentaion  with  absolute  positioning  using  the  TI  4100  and  the  GPS  system 
is  necessary  to  ascertain  the  accuracy  and  precision  with  which  solutions  can  be  made. 
Further  tests  should  be  made  to  compare  results  obtained  using  varying  collection  span 
periods  and  observation  intervals.  Improvements  in  satellite  ephemerides,  satellite  clocks, 
and  satellite  geometry  that  will  be  available  when  Block  II  GPS  satellites  are  deployed  may 
allow  absolute  positioning  with  sub-meter  accuracies. 
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